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We construct analytic formulas that represent the coupling-constant-averaged pair distribution 
function ~g xc (r s , C> fcpu) of a three-dimensional non-relativistic ground-state electron gas constrained 
to a uniform density with density parameter r a = (97r/4) 1 ^ 3 /fcF and relative spin polarization £ over 
the whole range < r s < oo and — 1 < £ < 1, with energetically-unimportant long range (u —* oo) 
oscillations averaged out. The pair distribution function g xc at the physical coupling constant is 
then given by differentiation with respect to r s . Our formulas are constructed using only known 
theoretical constraints plus the correlation energy e c (r s ,(), and accurately reproduce the g xc of the 
Quantum Monte Carlo method and of the fluctuation-dissipation theorem with the Richardson- 
Ashcroft dynamical local-field factor. Our g xc is correct even in the high-density (r s — > 0) and 
low-density (r s — > oo) limits. When the spin resolution of e c into ||, J.J., and |J. contributions is 
known, as it is in the high- and low-density limits, our formulas also yield the spin resolution of 
g xc . Because of these features, our formulas may be useful for the construction of density function- 
als for non-uniform systems. We also analyze the kinetic energy of correlation into contributions 
from density fluctuations of various wavevectors. The exchange and long-range correlation parts 
of our g xc (r s ,£,kFu) — 1 are analytically Fourier-transformable, so that the static structure factor 
S xc (r 3 ,C, k/kp) is easily evaluated. 



I. INTRODUCTION, DEFINITIONS, AND 
OUTLINE 

The exchange-correlation pair-distribution function 
g xc {v, r') of an iV-electron system is defined as 



9xc(r, r') 



N(N - 1) 
n(r)n(r') 



^ J |*(r, r',r 3 ...r J v)| 2 dr 3 ...dr A ,, 
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where n(r) is the electron density and '5 is the 
many-body wavefunction. Its coupling-constant average 
~g xc (r,r') is equal (in the Hartree units used throughout) 
to 



9xc( r , 



/ d\g* c (r, 
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where g xc {r, r') is the pair-distribution function when the 
electron-electron interaction is A/|r— r'| and the density is 
held fixed at the physical or A = 1 density. The coupling- 
constant averaged g xc plays a crucial role in density func- 
tional theory, .since it can account for the kinetic energy 
of correlation.!!! In fact, n(r') \g xc (v, r') — 1] is the density 
at r' of the exchange-correlation hole around an electron 
at r. 

In the uniform electron gas, n(r) — n and g xc (r.r') 
only depends on u = |r — r'|, and parametrically on 
the density parameter r s = (3/ 47m) 1 / 3 and on the spin- 
polarization ( — (N^rj- N±)/N. The coupling-constant 
average is in this casecl equivalent to an average over r s : 



9xc( r s,C,kFu) 



1 



9xc(r' s ,(,k F u)dr' s , (3) 



where kp = (97r/4) 1 / 3 /?' s is the Fermi wavevector. 
Clearly then 



d 

9xc(r s ,C,kFU) = ^—[r s g xc {r s ,C,y)] 
dr.. 



y—kp-u 



(4) 



and 



9xc( r ^C, k Fu) = gxc{Xr s ,(, k Fu). 



(5) 

The high-density (r s — > 0) limit is the weak-interaction 
limit in which the kinetic energy dominates. Relativis- 
tic effects are important for r s < 0.01. The low-density 
(r s — > oo) limit is the strong- interaction limit in which 
the Coulomb potential energy dominates. Fox r s > 100, 
the true ground-state density is not uniformJ3 but there 
is still a wavefunction that achieves the lowest energy of 
all those constrained to a given uniform density. 

The electron gas of uniform ^density is a paradigm 
of the density functional theoryEI for real, non-uniform 
electronic systems. The exchange-correlation energy of 
the uniform gas is the input to the local spin den- 
sity approximation, while the coupling-constant-averaged 
pair-distribution function is an input to the deriva- 
tion of gradient-corrected functionalsjffl to the construc- 
tion of the corresponding system-averaged exchange- 
correlation hole of a non-uniform density,u and to the 
implementation-.xiif. the fully-nonlocal weighted density 
approximationoQu We hope that our improved analytic 
model will be useful for these purposes, and also for 
the construction of new and more accurate functionals. 
In particular, the spin-resolved version of our model, 
when fully developed, could bring useful new informa- 
tion for the construction of functionals. Indeed, sim- 
ple hypotheses for the spin resolution have already been 
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used to construct several correlation functionals.i0 The 
uniform-gas a™ is also relevant to density matrix func- 
tional theory.ty 

The static structure factor S xc (r s ,£,k/kp) is the 
Fourier transform 

S xc (r s X,k/k F ) = 1 + — / {g xc (r s ,(,k F u) - 1] x 

, sm ku , , , , 

(kfu) — ; d(kpu), (6) 

ku 

and its coupling-constant average S xc is obtained by 
changing g xc into 7? xc in Eq. (||). Usually (fee and conse- 
quently <7 XC , S^c, and S xc are divided into exchange and 
correlation contributions: 

9xc{r s X, k F u) = g x (C, k F u) + g c (r s , C, fcpu), (7) 

where the exchange function g x is obtained by putting a 
Slater determinant of Kohn-Sham orbitals (or of Hartree- 
Fock orbitals) into Eq. ([!]). For a uniform electron gas, 
both Kohn-Sham and Hartree-Fock orbitals are plane 
waves, and g x is a simple function of kpu. The exchange- 
only pair-distribution function does not depend explicitly 
on r s , so that ~g x = g x : the explicit dependence on r a only 
appears when Coulomb repulsion is taken into account in 
the wavefuction. 

Both g x and g c have long-range oscillations. At high 
densities, these are Friedel oscillations; at low densi- 
ties, they represent the incipience of Wigner-crystal or- 
der within the liquid phase of uniform density. These 
oscillations are energetically unimportant in the follow- 
ing sensen A model which omits them but is constrained 
to have the same energy integral can correctly describe 
the short-range correlation while averaging out the oscil- 
lations of the long-range correlation. The energetic unim- 
portance of the oscillations is probably a consequence of 
the long-range and "softness'ppi the Coulomb interaction. 

Available analytic modelsoO of g c r&od, ~g c for the uni- 
form electron gas break down at hightijEJ (r s < 0.1) and 
low (r s > 10) densities. In this paper, we present a new 
model for the nonoscillatory part of g c (and hence g c ) 
which fulfills most of the known exact properties and is 
valid over the whole (0 < r s < oo) density range and for 
all spin polarizations Q. Our model is built up by inter- 
polatingbetween the short-range part recently computed 
in Ref. [lq and the long-range nonoscillatory part which, 
is exactly given by the random-phase approximation!!^ 
(RPA). Exact small- it and large- u expansions are recov- 
ered up to bishpr orders with respect to currently avail- 
able models.Bil3 All the parameters which appear in our 
interpolation scheme are fixed by exact conditions. We 
also build up a new nonoscillatory exchange g x which 
fulfills exact short-range and long-range properties up to 
the same order as our g c does. 

The paper is organized as follows. In Sec. |J, we ns t 
the known exact properties of g xc and g xc , and the ma- 
jor limitations of the models of Refs. || and [l2]. We then 



and for correlation (Sec. [V). In Sec. |v|, we discuss our 
results for exchange and correlation over the whole den- 
sity range. At metallic densities, we compare our analytic 
modelpwith the available Quantum Monte Carlo (QMC) 
dataj3ll3 finding fair agreement (Fig. 0). We also com- 



puted g c corresponding to t 
of Richardson and Ashcroft 



; dynamic local-field factors 
(RA), in order to see bet- 



ter how our model averages out the long-range oscilla- 
tions (currently not available from QMC). In this way, 
we are also able to show the effect of a dynamic local- 
field factor on the long-range oscillations, by comparing 
the RA result with the RPA (corresponding to zero local- 
field factor) long-range g c (Fig. ||). At high density, we 
find that our rnocjpl is in very good agreement with ex- 
act calculationsE3Ej (Fig. |J), and at low density it docs 
not break down and shows the expected (, dependence 
(Fig. |]). We also compare (Fig. ||) our model with pre- 
vious models|rta and discuss the qualitative effects of 
correlation (Fig. [?]). In Sec. VI, we discuss how to ex- 
tend our scheme to the spin-resolved (ft, II and T!) 
pair-distribution functions. The wavevector analysis of 
the kinetic energy of correl ation correspo nding to our S c 
and S c is presented in Sec. VII. Section VIII is devoted 



to conclusions and perspectives. 



II. EXACT PROPERTIES, AND LIMITATIONS 
OF PREVIOUS MODELS 



We list below most of the known exact properties of 
g xc and g xc for the 3D uniform electron gas. Equa- 
tion ([!]) implies the positivity constraint g xc > and 
the particle-conservation sum rule, which can be divided 
into exchange and correlation, 



du4iru 2 n (g x - 1) = -1 (8) 
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du Ami ng c = du Attu n g c 
'o 



(9) 



With the Coulomb interaction 1/u, the exchange function 
g x , the correlation function g c , and its coupling-constant 
averaged g c integrate to the exchange energy e x , to the 
potential energy of correlation v c , and to the correlation 
energy e c respectively, 

1 f 00 1 

- / duiTTu 2 -n(g x - 1) = e x (r s ,Q, (10) 

2 Jo u 



1 



1 



du4nu -ng c = v c (r s ,() (11) 



present our nonoscillatory model for exchange (Sec. Ill) 



1 f°° 1 

- / duATTu 2 -ng c = e c (r s ,Q- (12) 

2 Jo u 

For further discussion of the exchange hole density n {g x — 
1) surrounding an electron, the correlation hole den- 
sity ng c , and the generalization of Eqs. (||)-(|lj) to non- 
uniform densities, see Refs. ^| and |2^. 
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The short-range behavior of g xc is determined by the 
l/u Coulomb repulsion, which gives rise to the cusp 
conditioned 



dg x 



du 



u=0 



9x 



u=0 



(13) 



The function g xc satisfies a modified cusp conditiorU'El 
which can be derived from Eqs. (|^) and ([l3|). A quite 
accurate estimate of the r s and £ dependence of the 
short-range expansion coefficients of g xc and g xc has 
been recently obtained by solving a scattering problem 
in a screened Coulomb potential which describes the ef- 
fective electron-electron interaction, in a uniform elec- 
tron gas - the extended solutionis of the Overhauser 
model.tj (Classical electrons at zero temperature would 
have g xc \ u =o = 0, but nonzero values have a nondivergent 
potential-energy cost according to Eq. ([l|) and for quan- 
tum mechanical electrons lower the kinetic energy asso- 
ciated with the swerving motion needed to keep two elec- 
trons from colliding. Thus the right-hand side of Eq. ( |l3| ) 
is nonzero, except in the low-density limit. It is similarly 
nonzero for ia-,gas of classical electrons at an elevated 
temperatures) 

The long-range part of the nonoscillatory g xc corre- 
sponds to the small-fc behavior of the static structure 
factor, which is determined by the plasmon contribution, 
proportional to k 2 , and by the single-pair and multipair 
quasiparticle-quasihole excitaticp,CQntributions, propor- 
tional to k 5 and k 4 respectivelyESEj 



S xc {r s ,(,k -> 0) = 



k 2 



2ujp(r s ) 



0(k A ), 



(14) 



should hold for the ( ^ gas. Notice that the can- 
cellation of the-Jc 3 terms is obtained from beyond-RPA 
considerations. EZl 

Armed with these exact constraints, we can discuss the 
strengths and weaknesses of previous analytic . models, 
which unlike our present model break downtja outside 
the metallic density range 1< r s < 10. 

The Perdew-Wang modefl was largely based on first 
principles, plus limited fitting to Quantum Monte Carlo 
data. This model introduced the high-density limit of the 
RPA as the long-range component of g xc . But that limit 
was modelled crudely, leading to violation of the particle- 
conservation sum rule (and thus to failure for r s < 0.1). 
The model did not incorporate the plasmon sum rule, 
and produced an incorrect it -5 nonoscillatory long-range 
limit for g xc . The positivity constraint was violated at 
low densities, a problem evaded by switching over to a 
different analytic form for r s > 10. In this model, the 
spin resolution of g xc , even in its revised formjiil is less 
reliable than the total g xc . _ 

The model of Gori-Giorgi, Sacchetti, and Bacheletli£l 
was based upon extensive fitting to spin-resolved Quan- 
tum Monte Carlo data for £ = 0, and did not address 
nonzero £. Their model for g XCl unlike that of Perdew and 
Wang, was analytically Fourier-transformable to S xc . It 
incorporated the particle-conservation and plasmon sum 
rules, and the correct u~ s long-range limit for g xc , but did 
not build in the important high-density limit of the RPA 
for large u, leading to failure for r s <C 0.8. Moreover, 
small-u errors otthe Monte Carlo data were transferred 
into the modelEj 



where uj„(r s ) — y3/rf is the plasma frequency. Equa- 
tion (ph is called the plasmon swm rule. There is no 
k 3 term in the small-fc expansior£3 of S xc . Since, when 
k — > 0, the exchange-only static structure factor S x is 
equal to 



s.(c,*-o) = ! (i+o 2/3 + (i -o 2/3 



k 



k :i 



16k 3 / 



(15) 



there must be a linear term and a cubic term in the small- 
k expansion of the correlation static structure factor S c 
which cancel with the exchange. In real space, these 
terms correspond to long-range tails cx u~ A and oc u~ 6 
respectively.BES The nonoscillatory exchange-corpeJation 
pair-distribution function has a long-range taill3'c3 cx 
w -8 . As for more general densities, the exchange- 
correlation hole is more localized around its electron than 
the exchange hole (and thus better described by local 
or semi- local approximations for non- uniform densities). 
The high-density limit of thep^andom-phase approxima- 
tion (RPA) exactly describes^ the nonoscillatory long- 
range part of g xc , recovering Eq. ( pr|) through order k 2 . 
The absence of the A; 3 term in the small-A: expansion of 
S xc was demonstrated for the G-^j= gas by using ex- 
act frequency-moment sum rules.EH The same arguments 



III. NONOSCILLATORY EXCHANGE HOLE 

We present here our nonoscillatory model for the ex- 
change hole. This new model satisfies exact short-range 
and long-range conditions up to the same order as our 
correlation-hole model (Sec. IV) does. 

The exact exchange-only pair-distribution function for 
the uniform gas is 

g x (C,k F u) = l + i{(l + C)V[(l + C) 1/3 ^ U ] 

+(W)M(i-c) 1/3 fcH}, (16) 



where 



J{y) 



9 / sin y — y cos y 
2 I ? 



(17) 



Our nonoscillatory (J(y)) is parametrized as 



(J(y)) 



-9 



1 - e- A * y2 (l + A x y 2 + 



x y 



3! 



3-^ 2 (B x + C x y 2 + E x y 4 + F x y e 



(18) 



This model is similar in spirit, but not in detail, to those 
of Refs. U and ^8[ The first term of Eq. (|l^) achieves the 
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correct average long-range behavior — \y~ A as y — » oo, 
and is damped out at small y by the first square bracket 
which varies from y 8 as y — > to 1 as y — > oo. The 
second term then builds in the correct small-?/ behavior. 
The Gaussians smoothly blend the two terms, but are 
not motivated by any physical model. The analytic forms 
and linear parameters in Eq. ( |is| ) are convenient for con- 
straint satisfaction. The separation into long-range and 
short-range parts, although somewhat arbitrary, could 
be useful for the construction of new density functionals. 
The spherical Fourier transform of (J(y)), 



J{k) 



ky 



dy, 



(19) 



is also analytic and is reported in Appendix |a|. The large- 
ly expansion of Eq. (fill) is 



(Ay 



)) = -h- i + o(e-y ), 



(20) 



while the nonoscillatory average of the exact J(y) also 
contains a — f y -6 term (and no other long-range term). 
Such a term was included in the models of Refs. |^ and |2^, 
but with a coefficient wrong in both sign and magnitude. 
As explained in Sec. [n], the exact nonoscillatory corre- 
lation hole has long-range terms »~.f and y~ 6 which ex- 
actly cancel with the exchange ,OE3n3 so that the exact 
norioRcillatory exchange-correlation hole has a long-range 
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oc u which is purely correlation. However, 



tai 

as detailed in Sec. IV A, our nonoscillatory correlation- 
hole model is built without a u~ 6 long-range term, since 
this choice preserves a simple and useful scaling. We 
have thus also set the y~ 6 term to zero in our nonoscilla- 
tory exchange-hole model, in order to have an exchange- 
correlation hole with the exact u~ s long-range behavior. 

The six parameters A x through F x are fixed by requir- 
ing that (i) the particle-conservation sum rule is fulfilled, 
(ii) our g x gives zero contribution to the plasmon sum 
rule, (iii) our g x recovers the exact exchange energy, (iv) 
our g x is exact at u = in obedience to the Pauli prin- 
ciple in real space (two electrons of parallel spin cannot 
come together, since the antisymmetry of the wavcf unc- 
tion makes this probability vanish), (v) our g x has the 
exact second derivative at u = 0, and (vi) the informa- 
tion entropy S[— J(y)], 



poo 

S[-J(y)} = / dy4ny 2 J(y)ln[- 
Jo 



J{y% (2i) 



is maximized£3'E£l S of Eq. (|2^) is not a thermodynamic 
entropy but a mathematical one whose maximization en- 
sures that the analytic J(y) has no structure beyond 
that imposed by the exact constraints used to construct 
it. The parameter values are A x = 0.77, B x — —0.5, 
C x = -0.08016859, D x = 0.3603372, E x = 0.009289483, 
and F x = -0.0001814552. 

Our nonoscillatory model g x is compared with the ex- 
act exchange at £ = and £ = 1 in the upper panel of 
Fig. |IJ In the first panel of Fig. [|, the exchange hole 




2 2.5 3 3.5 4 

u/r s 

FIG. 1: Upper panel: our nonsoscillatory model for ex- 
change in the uniform electron gas is compared with the exact 
Hartree-Fock curve. Note that g x is the r s — » limit of g xc . 
Lower panel: low-density limit of our analytic model for the 
exchange-correlation pair-distribution function of the uniform 
gas. In this limit, the model g xc is almost exactly independent 
of the relative spin polarization (. 



g x — 1 is multiplied by {u/r s ) A in order to show how our 
model (solid line) averages out the oscillations of the ex- 
act exchange hole (dashed line). 



IV. NONOSCILLATORY CORRELATION HOLE 

Following Perdew and Wang^ we write the nonoscil- 
latory part of the correlation hole as the sum of a 
long-range part and a short-range part, somewhat as in 
Eq. M): 



(g c (r s X,k F u)} 



3 r s fi(v) 
k (kFu) 2 



[l-e"^ 2 ( 



1 + dx 2 



c n x n ~\ (22) 



where n = (4/3^)(9^/4) 1/3 , = [(l+C) 2/3 + (l-C) 2/3 ]A 
x = kpu/4>, and v = (j)K^/r^kFU. The six linear param- 
eters c n depend on both r s and Ci while the nonlinear 
parameter d only depends on £. 

The first term in the r.h.s. of Eq. (|22| ) is the long-range 
part of our g : the function fi(v) is a new parametriza- 
tion (sep-Sec. IV A ) of the RPA limit found by Wang and 
Perdewlij and displayed in Fig. 2 of Ref. ^. We multi- 
plied f 1 (v)/(kpu) 2 by a cutoff function which cancels its 
small-?i contributions, so that the long-range part of our 
g~ c vanishes through order u A and does not interfere with 
the short-range part. 
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FIG. 2: The function f{z, 0) given in Ref. [TJ]. The exact cal- 
culation (RPA) is compared with the present parametrization 
and with the one of Perdew and WangB (PW92). 



For modeling the short-range part, corresponding to 
the last term in the r.h.s. of Eq. (|22|), we use our im-. 
cent results obtained by solving the Overhauser model,t£l 
which allow us to fix the r s and C dependence of the lin- 
ear parameters c\, ci and C3 (Sec. IV B| ). We then use 
the remaining three linear parameters, C4, C5 and cq, to 
fulfill the particle-conservation sum rule and the plas- 
mon sum rule, and to recover the "exact" correlation en- 



ergy (Sec. IV C). Finally, the nonlinear parameter g?(£), 
which determines the "mixing" of long-range and short- 
range contributions, is fixed by impo sing the positivity 
constraint on g xc when r s — > 00 (Sec. IV D). 



A. Long-range part 

As discussed in Refs. |l6|and |[ the long-range (u — + 00) 
part of the nonoscillatory correlation hole can be ob- 
tained from the random-phase approximation by com- 
puting its r s — > limit. One finds 



n(g c (r s ,C, k F u)) -> ^ 3 (# 



Anv 2 



(23) 



where k s is the Thomas-Fermi screening wave vector, 
k s = Ky/r^kp. The function fi(v) is the spherical Fourier 
transform of the function f(z, 0) given by Eqs. (29), (34) 
and (36) of Ref. ^ 



f 1 (v)=2v 2 dzz 2 f(z,0) 



sin(wz) 



(24) 



where z — k/(f>k s is the proper scaled variable in recip- 
rocal space. The small- and large- z expansion of f(z,0) 
is 



f(z-> 0,0) = -- 
/(z->oo,0) = 



+ ^z 2 + 0(z 3 ) (25) 
0(z- 2 ). (26) 



2(l-ln2) „-l 

_9 * 



Equation (|26|) gives the high-density limit of the corre- 
sponding correlation energy, 



e c (r s - 0, C) = t^HO 3 lnr s + 0(r°), 




FIG. 3: Coulomb correlation contribution g c to the pair- 
distribution function g xc for the uniform electron gas for the 
paramagnetic (( = 0) and ferromagnetic (£ = 1) state. Our 
new analytic model is compared with the Diffusion-Quantum 
Monte Carlo results of Orjfeiz, Harris, and BalloneE-3 (OHB), 
and of Ceperley and Alderu (CA). The pair-correlation func- 
tion correspondingnto the local-field-factor model of Richard- 
son and Ashcroftcj (RA) is also shown. In the two bottom 
panels, the low-density limit of our g c is reported. 



which is exact at £ = and 1, but is slightly differ- 
ent from the exact result for < Q < 1 (see Refs. ^ 
and [l6| for further details). The small- z expansion of 
f(z,0), Eq. (|25|), fulfills the particle-conservation sum 
rule [f(z = 0,0) = 0], contains a linear term which can- 
cels with the exchange (and corresponds to a long-range 
tail oc w~ 4 in real space, see Sec. ||), and fulfills the plas- 



mon sum rule [exact z 2 coefficient, see Eq. (14)]. The 



term in Eq. (E5I), if it does not vanish, produces a u 6 
contribution to the correlation hole at large u. 

As said in Sees. |n| and III, the long-range (u — > 
00) nonoscillatory behavior of the exact exchange 
hole coniaiii£L ?'~ 4 and u~ 6 contributions which are 
cancellccoEj'EZl by similar contributions to the exact cor- 
relation hole. When we use the high-density limit of 



(27) 



Eq. (23) for the long-range part of the correlation hole, 



we automatically achieve cancellation of the u terms. 
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But to cancel the u 6 terms in g x — 1, we would have to 
replacejf^u)/!- 2 in Eq. @ by Ji{v)/v 2 + r s cj)h{r s , C, v), 
where f 1 (v)/v 2 has no v contribution and h is propor- 
tional to v~ e with no r s or £ dependence at large v. The 
extra term r s cf>h vanishes in the high-density limit for a 
given v, and is unknown. Since we want to keep for our g c 
the simple form of Eq. (22), but we also want to have the 
correct long-range behavior (oc u~ 8 ) for ~g xc , we decided 
simply to set the u~ e terms to zero in both our exchange 
(Sec. HI) and correlation- hole models. Figures || and^j 
do not suggest that this choice introduces any significant 
error into our models for the separate exchange and cor- 
relation holes. 

We thus parametrize f 1 (v) as follows 



fi(v) 



a + b 2 v + a\v 2 + a 2 v A + a 3 v 6 
{v 2 + b 2 ) 4 



(28) 



With respect to the parametrization given by Perdew and 
Wang,H our Eq. ( p8f ) has the advantage that it is analyt- 
ically Fourier-transformable (see Appendix |b]), so that 
the particle-conservation sum rule and the plasmon sum 
rule can be easily imposed. (They are not fulfilled by the 
Perdew and Wanga parametrization). After imposing on 
our f 1 (v) all the exact properties plus the vanishing of 
the z 3 term in Eq. (^H|), we are left with one free parareu 
eter, b, which is fixed by a best fit to our RPA data.Ej 
All the parameter values are reported in Appendix 
The function f(z , 0) corresponding to our parametriza- 
tion [see Eq. (p3l[)] is compared in Fig. || with the RPA 
result|-and with the Fourier transform of the Perdew and 
WangB (PW92) J^v). 



B. Short-range part 



Our g c has the small-u expansion 



k,FU , , . ( Ufu 
\9c) = ci + c 2 — — + (-cid + c 3 ) 



0(u 3 



(29) 



In order to recover the short-raage behavior obtained by 
solving the Overhauser modeljlj we require 



(l-C 2 ) 



(30) 



C2 = ^(^) 1/3I ¥ !l f(1+C)1/3 t (1 - C)1/3l ^ i H i ) (3D 

C3 = ^ 2 [(^) 2/3 ^(^.c)- (i+o8/ ;r' c)8/3 ] + 

+cid, (32) 



where rj* = 2r,/[(l + C) 1/3 + (1 - C) 1/3 ], and aj 1 , a\ l 
and a 2 (r s ,^) are given by Eqs. (36), (37) and (46) of 
Ref. [l^. In this, way, the modified cusp condition is ex- 
actly satisfiedJlj 
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FIG. 4: Upper panel: long-range part of the exchange hole. 
Our nonsoscillatory model is compared with the exact ex- 
change. Second and third panel: long-range part of the cor- 
relation hole. Our nonoscillatory model is p«mpared with g c 
obtained from the Richardson and Ashcrofta (RA) local-field 
factor. In the lowest panel the random-phase-approximation 
(RPA) result for r s = 10 is also shown. All curves are for the 
C = gas. 



C. Sum rules 

We want our correlation hole to satisfy the particle- 
conservation sum rule and the plasmon sum rule, 
and to recover the "exact" correlation energy. Our 
new parametrization of the function fi(v) satisfies the 
particle-conservation sum rule, and recovers the exact 
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plasmon coefficient and the lnr s term of the resulting 
correlation energy. Thus, we only have to require that 
the remaining part of our ~g c gives zero contribution to 
(i) the particle-conservation sum rule and (ii) the plas- 
mon sum rule, and (iii) recovers the correlation energy 
beyond the lnr s term. In this way we have three linear 
equations for the three parameters C4, C5 and Cq: 



6 />oo 

Va„ / e- t2 t n+1 dt = AS{a) 

n=l Jo 
6 p 00 

Y\cn / e- t2 t n+3 dt = AP{a) 
n=l J ° 

6 r oo 

n=l 



(33) 
(34) 
(35) 



n 



where c„ = c n /d ^ , t — y/dkpu/cf), A = (f>r s d/K, a — 
</) 2 K(r s /d)^ 2 , and 

poo 

S{a) = / Jiiat^ 2 (1 + t 2 + i< 4 ) dt (36) 



P(a) 



f 1 {pd)e- t t 2 (l + t 2 + i< 4 ) (37) 



i?(a) 



/i(at) 



l-e _t (l + t 2 + it 4 ) (38) 



2/3 



The functions 5 f (a), P(a) and i?(a) can be obtained an- 
alytically and are reported in Appendix |^. The param- 
eters C4, C5 and cq are then equal to 

5 4 = {100V^(3tt - 8)ci + (690tt - 2048)5 2 + V¥(225?r 
-672)5 3 + (8192 - 2100n) AS (a) + AP(a){600n - 
2048) + 960V7r[AP(a) - E]}/[4(512 - 165vr)] (40) 

c 5 = 2{(30tt -128)5i -8V^c 2 + (397T- 128)5 3 - 
lUy^AS{a) + 16y/rAP(a) - 256[AR(a) 
-£]}/(512-165tt) (41) 

ce = {^(18071 - 624)5i + (150tt - 512)2 2 + V^(135tt 
-432)5 3 + (3072 - 1260tt) AS '(a) + AP(a){360n - 
1024) - 480y/^[AR(a) - P]}/[6(165tt - 512)]. (42) 

D. Positivity constraint in the low-density limit 

The nonlinear parameter d can be fixed by imposing 
the condition that g xc remains positive when r s — > 00. 
The short-range behavior imposed on our g c ensures that 
the small- w expansion of the corresponding ~g xc has coeffi- 
cients which are always > through order u 2 , and which 
become zero in the low-density or strongly-correlated 
limit. We have checked that, if we want to have a pos- 
itive g xc for all densities, we only need to require that 
also the u 3 coefficient (equal to C4 — dc2) becomes when 




FIG. 5: Coulomb correlation contribution to the pair- 
distribution function for the uniform electron gas for the para- 
magnetic (£ = 0) and ferromagnetic (£ = 1) state in the high 
density [r 3 — * 0) limit. The result from our analytic model is 
compared with the exact calculation of Refs. |13|Jl9| . 



r s —> 
pairs 



fording to the cusp condition for parallel-spin 
- : ' ^ We thus have an equation for d((): 



lim c 4 (r s ,()-d(()c 2 (r s ,() = 0. 



(43) 



Equation ( f43| ) is rather complicated since C4 also depends 
nonlinearly on d. However, it can be sobzed numerically 
for each £, and, when the Perdew-WangLj parametriza- 
tion of the correlation energy is used in Eq. (|39|), the 
result is very well fitted by 



d(C) = d(o) [(i + C) 2/3 + (i-C) 2/3 
with d(Q) = 0.131707. 



(44) 



V. RESULTS FOR THE 
EXCHANGE-CORRELATION HOLE 

In the next three subsections we present and discuss 
our results for the nonoscillatory g x , g c and g xc in the 
whole (0 < r s < 00) density range. We have used the 
correlation energy e c as parametrized by Perdew and 
Wang,E2l which was built with the Quantum Monte Carlo 
data of Ref. [| as an input. It is however straightforward 
to build into our equationiLan ab initio e c for the 3D uni- 
form gas when available£j showing that the exact con- 
straints suffice to determine g xc without the need for any 
"numerical experiment" . 



A. Metallic densities 

In the six upper panels of Fig. || we compare our an- 
alytic g c with the Quantum Monte Carlo (QMC) data 
of Ceperijey and Alden3 (CA) and of Ortiz, Harris and 
BallonetZl (OHB) for r s = 2, 5 and 10, and for C = 
(left) and £ = 1 (right). In the £ = case, we also report 
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g c as obtained by the dynamic local-field-factor model 
of Richardson and AshcroftJia (RA) via the fluctuation- 
dissipation theorem (as in Ref. [35]) . The RA model yields 
very accurate correlation energies e c (r s , £ = 0),c3 and we 
find that the RA g c is in very good agreement with QMC 
data except at small u. The limit u — ► is not correctly 
included in the RA parametrization of the local-field fac- 
tor, which violates the Pauli principle in real space. 

We see that our model is in fair agreement with QMC 
data for the paramagnetic gas. In the ferromagnetic case, 
where the pair-correlation function shows stronger oscil- 
lations even at intermediate densities, the agreement is 
less satisfactory (as in the model of Ref. [2]). This is not 
surprising, since our model does not take into account the 
energetically unimportant oscillations: it only includes 
the minimum number of oscillations needed to fulfill the 
sum rules. This is evident in the second and third panel of 
Fig. ^, where g c is multiplied by (u/r s ) 4 . In this way, the 
long-range oscillations are amplified and become clearly 
visible even at metallic densities: we can thus compare 
our model (solid line) with the RA result (dashed line). 
This is done at r s = 2 and 10. The many exact prop- 
erties imposed on the RA local-field factor and the first 
three left panels of Fig. || suggest that the long-range 
part of the RA g c is very reliable and that the oscilla- 
tions are probably accurately described. One clearly sees 
in Fig. |] how our model follows the first oscillation and 
averages out the others. In the lowest panel of Fig. ||, 
the long-range oscillations of the random-phase approxi- 
mation (RPA) g c at r s = 10 are also shown. At large r s , 
the RPA oscillations of g c tend to cancel the ones of g x 
(first panel), while the effect of a dynamic local-field fac- 
tor clearly inverts this tendency: the oscillations of the 
RA g c (second and third panel) are almost in phase with 
the oscillations of g x . We interpret this to mean that the 
RA g xc of the low-density uniform electron gas is build- 
ing up an incipient Wigner-crystal-like order of the other 
electrons around a given electron. 



B. High density 
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FIG. 6: Comparison of the present work with the models 
of Ref. | (PW92) and of Ref. |l| (GSB) at high densities 
(first panel), metallic densities (second panel) and in the low- 
density regime (third panel). In the r s — 100 case the original 
PW92 curve has been divided by 10, and the low-density form 
proposed in the Appendix of Ref. (PW92-App) is also re- 
ported. All curves are for the paramagnetic (£ = 0) gas. 



In the high-density limit, g c qi 27} \. goes to zero, so that 
9xc — * 9x- It has been shownoEjO that in the r s — > 
limit g c /r s remains finite and goes to a well defined ianc- 
tion of u/r s , which has been computed exactlyJiSlij In 
Fig. [5] we compare this exact calculation (dashed line) 
with our model (solid line), computed at r s = 10~ 5 , for 
( = and £ = 1. We see that (i) our model does not 
break down as r s — > 0, and (ii) there iSi-£air agreement 
with the exact result. Previous modelsota for g c usu- 
ally break down at r s ~ 0.1. Feature (i) is due to the 
new parametrization of fi(v) which exactly fulfills the 
particle-conservation sum rule, while feature (ii) is due 
to the short-range behavior taken from Ref. 15 which 



includes the exact high-density limit of the short-range 
coefficients. 



C. Low density 

In the low-density or strongly-correlated limit, we ex- 
pect that g xc (equal to g xc in this case) does not depend 
on since in this limit the Pauli principle in real space 
becomes irrelevant with respect to the Coulomb repul- 
sion. In the lower panel of Fig. [j] we report our model 
at r s = 10 5 for three different values of the spin po- 
larization £. We see that the C dependence of our low- 
density g xc is indeed H#ty weak, and that, unlike previ- 
ous parametrizations,ali3 our model never gives rise to 
an unphysical negative pair-distribution function. Fig- 
ure [l] also offers a view on the same scale of the extreme 
high-density limit of g xc (equal to the exchange-only pair- 
distribution function, first panel) and of the extreme low- 
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E. Features of the "correlation factor" 
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FIG. 7: The correlation factor K xc defined in Eq. (]4^) for the 
paramagnetic (upper panel) and ferromagnetic (lower panel) 
uniform electron gas. 



density limit (second panel). We see how the ( depen- 
dence of g xc , which is very strong in the r s — > limit, is 
cancelled by correlation in the r s — > oo limit. The low- 
density limit of our g c = g xc — g x is reported in the two 
lowest panels of Fig. || for £ = or £ = 1. 



D. Comparison with previous analytic models 

In Fig. |^ the present model is compared with the 
parametrizations of Perdew and Watigcl (PW92) and of 
Gori-Giorgi, Sacchetti and Bacheletli3 (GSB). In the first 
panel, we see that in the high-density cesimc (r s = 0.01) 
the PW92 model starts to break down,ESEj and that the 
GSB parametrization is completely unable to describe 
such high densities. (This is due to the wrong r s — > be- 
havior of the GSB on-top pair density.) At r s = 2, well 
inside the metallic regime, we see (second panel) that the 
present work is very close to the PW92 model and slightly 
deviates from the GSB curve at u/r s < 1. Finally, in the 
third panel we show the total pair-distribution function 
g xc at r s = 100: the PW92 model in its original form 
completely blows up, while the GSB model becomes neg- 
ative at u/r s < 1 but is still "reasonable". The low- 
density form proposed in the Appendix of Ref. ^ (PW92- 
App) is also reported: it corresponds to an exchange- 
correlation hole narrower than the present one. 



To better see the effects of correlation, we define a 
"correlation factor" 



K xc (r a X, kFu) 



9xc 



1 



9x 



1 



1 



9c 



9x 



1 



(45) 



which morphs the exchange hole into the exchange- 
correlation hole, and is displayed in Fig. ^. We must 
of course use non-oscillatory models here, since the ex- 
act g x — 1 has nodes which would create singularities 
in Eq. (^5|). Figure (?] shows that K xc — > K x = 1 in 
the r s — > limit. For typical valence-electron densi- 
ties, we see that correlation enhances or deepens the hole 
(K xc > 1) around an electron for u/r s < 1.5, while it 
screens out the long-range part of the hole. Because of 
the exact cancellation of the u -4 and u~ 6 long-range 
terms between g c and g x — 1, K xc at large u goes to 
like u -4 . For r s > 2, K xc can be negative in the 
range 1.5 < u/r s < 3, corresponding to a positive peak in 
g~ xc — 1. We can think of K xc (r s , £, kpu)/u as an effective, 
density-dependent screened electron-electron interaction 
whose exchange energy equals the exchange-correlation 
energy of the Coulomb interaction 1/u. . .. . 

The correlation factor has a possible applicationEllEJ 
to the modelling of exchange and correlation in systems 
of non-uniform density. First we note that the exchange- 
correlation energy is fully determined by the spherical 
average n xc (r, u) of the hole,S 

n xc (r,u) = J -^-n xc (r,r + u). (46) 

A possible "correlation factor model"0 for n xc (r,u) is 
n xc (r, u) = K xc (r, u)n x (r, u), (47) 

where ^(r, u) is the exact exchange hole. K xc for a non- 
uniform density could be constructed from Eq. ( (45| ) by 
inserting into Eq. ( p2| ) (or a simplification thereof) an in- 
dependent set of linear parameters c n chosen to satisfy 
exact constraints on n xc . The result would presumably 
be a model for exact exchange and approximate corre- 
lation compatible therewith. The screening of the long- 
range part of the exact exchange hole is essential for a 
proper description of molecules.E3 

VI. SPIN RESOLUTION 

We can define spin-resolved pair-distribution functions 
which describe spatial correlations between ||, J. J,, and 
t| electron pairs. Their normalization is such that the 
spin-averaged g xc of Eq. (|l]) is equal to 



i + C 



>Jxc 



i-C 



(48) 
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While the spin resolution of the exchange-only pair- 
distribution function g x is well known,0 the correlation 
part is much more delicate, and an accurate analytic rep- 
resentation is only available^ for ( = in the density 
range 0.8 < r s < 10. 



The model presented in Sec. |IV|can be used to build up 
spin-resolved correlation functions provided that the spin 
resolution of the input quantities is known. The input 
quantities are (i) the RPA long-range part, (ii) the short- 
range coefficients from the solution of the Overhauser 
model, and (iii) the correlation energy. Once these input 
quantities are known, in fact, one can build, say, g~1^, 
starting from the same Eq. (p2|), using the RPA fj, long- 
range part, and putting the f J, short-range coefficients 
§, and e rl into Eq. ©. Finally, the 



into Eqs. 

positivity constraint of g\} c in the low-de nsity limit can 



IV D 



be applied to find dfj.(C)j as done in Sec. 

The first point is thus to see whether the quantities 
(i)-(iii) are available in their spin- resolved contributions. 
The RPA long-range part is easily spin-resolved for the 
C = gas,Ejo while its spin resolution in the partially 
polarized gas is less trivial. The short-range coefficients 
from the Overhauser model are available as ||, J, J, and 
t| separate contributions.E3 The correlation energy rep- 
resents the major problem: at ( = it can be easily 
spin resolved in the high- and low- density limits, while 
at intermediate densities the best estimate is probably 
the one given in Ref. [l2| Almost nothing about the spin 
resolution of e c is known for the ( ^ gas, except in 
the extreme low-density limit, when the system becomes 
(■-independent. 

Here, we show results for g]} in three cases: the ex- 
treme low-density limit, the high-density limit of the 
paramagnetic gas, and the r s — 2, ( — case. The low- 
density limit must be treated first, since it is necessary 
to determine d|j.(C) through the positivity constraint on 
g].\, when r s — > oo. 

When £ = 0, the spin-resolution within RPA is very 
simple: up-up and up-down inter actions contribute the 
same amount to correlation rWi The long-range part 
of our gj} can thus be built using the function fi(v) 
of Eq. ( p8| ) with the same parameters of Appendix M. 
While the spin-averaged nonoscillatory long-range behav- 
ior computed within RPA is also exact beyond it at all 
densities, its sain resolution is exact beyond RPA only 
when r s — > 0.L3 We keep on using it even in the extreme 
low-density limit, since it is the only way to build up a 
spin-resolved g c starting from our model. As we shall 
see, the results obtained are reasonable, and justify our 
choice. When r s — > oo, we. expect the statistics to be 
energetically unimportant so that e x \ = e x \, — = 
e xc . We thus find e\} — e xc — — 0.892/r s , where the 
numerical coefficient corresponds to the Perdew-Wang 
parametrizatioro of e c . The positivity constraint on g x ^ 
gives 



%(C) = <%(0) (i + C) 2/3 + (l-O^-i 



\2/3 



(49) 




FIG. 8: Up-down pair-distribution function for the uniform 
electron gas in the low-density limit. 
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with d n (0) = 0.0885717. The results for g^ c are shown 



FIG. 9: Up-down Coulomb correlation contribution to the 
pair-distribution function for the paramagnetic (£ = 0) uni- 
form electron gas in the high density (r s — > 0) limit. The 
result from our analytic model is compared with the exact 
calculation of Ref. FL9l 



in Fig. ||, at r s = 10 5 , for £ = and £ = 1. 

For the high-density limit of the paramagnetic gas, all 
the spin-resolved input quantities are exactly known. It 
is thus the best case to test our model. When r s — > 0, 
the spin-resolution from RPA is exact also beyond it: the 
long-range part of gj* is in this case exactly described 
by Eq. ( |28| ) with the parametets-pf Appendix |^. The 
correlation energy, in this limit ,E3c3 is simply equal to the 
spin-averaged correlation energy of Eq. ( p7| ) with £ set to 
zero. The short-range || coefficients from Ref. [TBI include 
the exact spin-resolved high-density limit of the ( = 
gas. The so-obtained g]} is shown in Fig. || together with 
the exact calculation from Ref. [l^, which is, in this case, 
equal to the RPA result. We find very good agreement. 

At metallic densities, we used the spin-resolved e c for 
the C = gas from Ref. [l2|, and the RPA spin resolution 
for the long-range part. In Fig. [Io[ we report our results 
for g£I anc i gV\ = 2g^ — at r s = 2, together with the 
QMC data of Ref. |l^, and with the values that we have 
obtained frop, the Richardson and Ashcroft (RA) local- 
field factors. lia We see that our result is reasonable, but 
does not accurately agree with the QMC data. In this 
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FIG. 10: Spin-resolved Coulomb correlation contribution to 
the pair-distribution functions for the paramagnetic uniform 
electron gas at density r a — 2. The present model is compared 
with the Quantum Monte Carlo (QMC) data from Ref. |l| 
and with the result obtained from the local-field factors of 
Richardson and Ashcroft£3 (RA). 



respect, the RA results are much better for u/r s > 0.7, 
while they blow up in the short-range part, since they do 
not satisfy the Pauli principle in real space. As said, the 
spin resolution is very delicate, so that an analytic model 
is very difficult to build up. The best analytic represen- 
tation of gl^ and gl' at metallic densities is probably the 
one of Ref. 12, which was built to interpolate the QMC 
data of Ref. |17| accurately. 



VII. WAVEVECTOR ANALYSIS OF THE 
KINETIC ENERGY OF CORRELATION 

Wavevector analysis^ is usually a study of the static 
structure factor of Eq. (^J). The wavevector analyses of 
the correlation energy e c and of the potential energy of 
correlation v c have often been reportedJ13c3'E3 while the 
kinetic energy of correlation t c is much less studied. We 
can decompose t c (r s , £) into contributions from different 
wavevectors of a density fluctuation, 



/>oo 

tc(r.,C) = | / dqq 2 T c (r s X,q) 1 
Jo 



(50) 



where q = k/kp. Since t c — e c — v c , the wavevector 
analysis of t c is just the difference between those for e c 
and v c - 

rn > „s 2fc F [S c (r s ,t,q) ~ S c (r s ,(,q)} 

%{r s , C) Q) = 1 5 ' 51 ) 

d7T q z 

The small-g limit of T c can be obtained by the plasmon 
sum rule, 



T c (r s ,C,9^0) = ^r7 3 / 2 + O(g 2 ), 



(52) 



and its leading term is independent of £, as expected 
from Eq. (|l4|). To write down the large-g limit of T c we 



need^p-epaiand S c and S c for large arguments. We know 

thateeBS 



S c (r s ,(,q^> oo) 



4 2g xc (r s ,(,u = 0) 



Snki 



+ 0(q- 6 ), 
(53) 



from which we can also obtain the large-g limit of S c 



where 



S c (r s ,(,q 



1/3 



(x ) = l i + 0{q- & ), 



(54) 



7=-=- cT -/ r' s g xc (r' s ,C,u = 0)dr' s . (55) 
37T \9n J r s J 

Through the cusp condition of Eq. ([l3"|), we see that the 
large-g limit of S c is determined by the coefficient of u/r s 
in the small-u expansion of ~g c , ai(r s ,£) [see Eqs. (35) 
and (37) of Ref. f|, related to c 2 (r s ,0 of Eq. @ by 
ai = (97r/4) 1 / 3 c 2 /0. We thus have 



T c (r s ,(,q^> oo) = 



2g xc (r s ,C,u = 0) 



)Qi(r s ,C) 



9tt 2 



(56) 

In Fig. |ll|, we report 7^ for the C — gas, for two differ- 
ent densities, r s = 2 and r s = 5. We clearly see that the 
small wavector contribution to t c comes from the kinetic 
energy of the long-wavelength zero-point plasmons, and 
that the decay of the plasmon contribution with increas- 
ing wavevector k is gradual. The corresponding result 
from the Richardson and Ashcroft local field factonla is 
also shown. The Richardon- Ashcroft model gives-a good 
description of plasmon dispersion and dampingo 

It is also interesting to compare T c with the decomposi- 
tion of the kinetic energy of correlation into contributions 
from different wavevectors of a quasi-electron. For £ = 0, 
we can write 



< c (r s ,C = 0) = 



dqq 2 n c (r s ,C = 0,q)(k F q) 2 . (57) 



Here n c is the correlation contribution to the momentum 
distribution, n c (q) = n(q) — no(q), and no is the Fermi 
step function. The leading term in the small-g expansion 
of kpq 2 n c is proportional to q 2 , and is thus rather differ- 
ent from the corresponding behavior of T c , Kq. (^2|). On 
the other hand, in the large-g limit we havea 



k 2 F q 2 n c (r s ,( = 0, q -> oo) 



9xc(r s ,C = 0,u = 0) 



9tt 2 



q 



(58) 



a behavior very similar to Eq. (|56j). This is not surpris- 
ing, since the large-g limits of both S c and n c are deter- 
mined by the downward-pointing kinks in the many body 
wavefunction [producing the cusp of Eq. ([r])] which oc- 
cur whenever two electrons of antiparallel spin come to- 
gether. In the r s — > limit, Eqs. ([56|) and ( |58| ) become 
equal. A study of the equations linking S c and n c from 
the point of view of density matrix functional theory is 
reported in Ref. Ell 
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FIG. 11: Wavevector analysis of the kinetic energy of corre- 
lation at r a — 2 and r s — 5 for the paramagnetic gas. The 
function T c {r s ,C,,k) is denned in Eq. (jHl]). The present work 
is compared with the result obtained from the Richardson and 
AshcrofO (RA) local-field factor. 



VIII. CONCLUSIONS AND FUTURE 
DIRECTIONS 



wavevector analyses of the kinetic energy of correlation in 
the uniform electron gas, finding them the same only in 
the limits of large wavevector and high density. We have 
also found that the decay of the plasmon contribution 
with increasing wavevector fc is gradual. 

In the future, it may be possible to construct the cor- 
relation energy e c (r s , £) and its spin resolution eJ2L(r s , Q 
directly by interpolation between known limits, lj with- 
out using any Monte Carlo or other data. This develop- 
ment would probably not give us a better e c (r s ,£) than 
we already have, but would provide the first spin reso- 
lution over the whole range of r s and £; it would also 
show that the known exact constraints are by themselves 
sufficient to determine g xc . The extended Overhauser 
modeluJ might be evaluated for ( different from zero, to 
test and refine the spin-scaling relations used in Ref. 
The extendedpOverhauser model can also be made more 
selfconsistent .c3 

A small FDRTRAN77 subroutine which numerically 
evaluates our g r [Eq. (|22|)1 can be downloaded a t 
tittp : //axtnt2 . phys . uniromal . it/PGG/elegas .html. 



The known exact constraints summarized in Sec. ||, 
plus the random phase approximation for long-range 
(u — > oo) correlation, the extended Overhauser modelMI 
for short-range (u — > 0) correlation, and the correla- 
tion energy e c (r S) 0, suffice to determine the pair dis- 
tribution function g xc (r s , kpu) of a uniform electron 
gas over the whole density range, including the high- 
density (r s — » 0) and low-density (r s — > oo) limits, apart 
from energetically-unimportant long-range oscillations. 
The analytic formulas we have so constructed for the 
coupling-constant-averaged g xc should be useful for fur- 
ther developments and applications of density-functional 
approximations for the exchange-correlation energy of a 
non-uniform density. 

For metallic densities (2 < r s < 10) with £ = or 
1, ourj,Cx|c is in good agreement with Quantum Monte 
CarlcEO In the same density range for £ = 0, it 
also agrees with the-@ xc we have calculated from the 
Richardson- Ashcroftlia dynamic local-field factor, except 
near u = where the Richardson- Ashcroft model was 
found to break down (although this model seems to 
describe the long-range oscillations correctly). In the 
r s — > limit for small k F u 1 our g xc agrees with the re- 
sults pfjaerturbation theory to zero-th (exchange) or first 
orderEjtJ in the electron-electron interaction. The static 
structure factor S xc is also modelled accurately, neglect- 
ing the non-analytic structure of the.exact S xc at k = 2kF 
arising from long-range oscillations r3 

Our formulas can also be used to spin-resolve g xc into 
|f, H, and t| components (Sec. |v| ), when the spin reso- 
lution of e c is known (as it is in the high- and low-density 
limits) . The additional information in the spin resolution 
might well be used to construct more accurate density 
functionals for the correlation energy. 

We have also examined two physically-different 
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APPENDIX A: NONOSCILLATORY EXCHANGE 
HOLE IN RECIPROCAL SPACE 

In reciprocal space, the exchange-only static structure 
factor is equal to 



S x (C,k/k F ) = 1 + — 

6TT 



(1 + QJ 



k 



Mi + C) 1/3 



(i-C)J 



MW) 1/3 



(Al) 



where J(k) is defined by Eq. (|19[). From our parametriza- 
tion of (J(y)) [Eq. @] we obtain 



J(k) 



32 C 



k 2 - 6A X \ y/n _ k 2 




D 



.3/2 



C x (6D x ~k 2 ) E x (60Dl -20D x k 2 + fc 4 ) 



4Dl /2 



F x (MQDl - A2QD 2 x k 2 + A2D x k i - fc 6 ) 



MD 



15/2 



(A2) 
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APPENDIX B: LONG-RANGE CORRELATION 
HOLE IN RECIPROCAL SPACE 



We obtain for S(a), P(a) and R(a): 



f(z,0) 



1 



The function f(z, 0) corresponding to our Eq. ( p8|) is S(a) — aol2(a) + (ao + aia 2 )Tf(a) + (|ao + a\a 2 + 

a 2 a 4 )2" 4 (a) + (^aia 2 + a 2 a 4 + a 3 a 6 )2" 4 (a) + 
(ia 2 a 4 + a 3 a 6 )2f (a) + ia 3 a 6 X 4 10 (a) + b 2 a x 
[Ji(a)+22(a) + ^j|(a)] (C3) 
3a 2 6 6 + 9a 3 6 8 )z 2 + (a 6 3 - ai b 5 + a 2 b r p( a ) = aoZf (or) + (ao + aia 2 )2" 4 (a) + (|a + ai« 2 + 



{ 



, -si- — [48a + (33a fc ™ 3ai6 3 
2zb 8 I 48 L 

-3a 2 6 5 - 15a 3 6 7 )z + (9a 6 2 - 3ai6 4 - 



-a 3 &V]} 



6 2 d 3 X(6,z) 
6ttz d(b 2 ) 3 ' 



where 



l(b,z) = - [e hz E x {bz) ~ e- bz E x {-bz)} , 
and, with x > 0, 



(Bl) 



(B2) 



#i(aO= / —dt 

J x t 

Ei(-x) = -Ei(x) = -PV (^J jd?j . 

Here PV means the Cauchy principal value integral.! 4 ! 

The parameter values which satisfy Eqs. (|2^) and (p6[). 
give rise to a zero coefficient for the z 3 term in the smallpZ. 
expansion of f(z,0), and accurately fit our RPA data,E£l 



no = 26 8 C 
66 3 

oi = ^ r k 2 Co6 3 + 786- 256v / 3l 

a 2 = 486 2 /vr 2 
a 3 = 12/tt 2 
36 4 

b 2 = — 



C 
6 



■ [96^3 - 366 - 6 3 C* 7T 2 
-2(l-ln2)/7r 2 



7. 



R(a) 



a 2 a 4 )Zf (a) + {^a ia 2 + a 2 a 4 + a 3 a 6 )2f (a) + 
(±a 2 a 4 + a 3 a 6 )ll"(a) + ±a 3 a 6 2" 4 12 (a) + b 2 a x 
[Z 3 (a)+Z 4 5 (a) + ±Z 4 7 (a)] (C4) 
a l^ 1 (a) - (ao + aia 2 )2" 4 (a) — (|ao + aia 2 + 
a 2 a 4 )2f (a) - ( ±aia 2 + a 2 a 4 + a 3 a 6 )2"f (a) - 
(^a 2 a 4 + a 3 a 6 )Xj(a) — 5a 3 a 6 X|(a) — 6 2 a x 
[X 4 °(a)+X 2 (a) + iX 4 (a)] + 



2ai + a 2 6 2 + 2a 3 6 4 5 7T 
1266 + 2 32 67' 



(C5) 



where ao, ai, a 2 , a 3 , 6 2 a nd 6 are given in Appendix 
The integrals of the kind (CI) can be written as 



2£(a) = ^(«/6)/6 2 ", 



o [(rx) 2 + l]' 



-dx, 



(C6) 
(C7) 



and starting from 



2(0 - ^ 1/r 



1 - erf I - 

r 



(C8) 
(C9) 



APPENDIX C: ANALYTIC EXPRESSIONS FOR 
THE FUNCTIONS 5(a), P(a) AND 72(a) 

The three functions of Eqs. (|36|), ( ^7| ) and (|38|), which 
enter our model for <7 C , are given by a linear combination 
of integrals of the kind 



o [(ax) 2 + 6 2 ]' 
1 - e-* 2 



o x[(ax) 2 +6 2 ] 



■dx 
— dx. 



(CI) 
(C2) 



can be computed by different iati on with respect to a and 
6. The integrals of the kind (C2) can be also obtained by 
differentiation with respect to 6 of 



r-l M e b2 '°< 2 E l {b 2 /a 2 ) + Hb 2 /a 2 )+ 1 

T i w = ^ . (CIO) 



where 7 = 0.5772156649. 
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